Generic Phase Diagram of Fermion Superfluids with Population Imbalance 
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It is shown by microscopic calculations for trapped imbalanced Fermi superfluids that the gap 
function has always sign changes, i.e., the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state like, up 
to a critical imbalance Pc^ beyond which normal state becomes stable, at temperature T = 0. A 
phase diagram is constructed in T vs P, where the BCS state without sign change is stable only 
at T 7^ 0. We reproduce the observed bimodality in the density profile to identify its origin and 
evaluate Pc as functions of T and the coupling strength. These dependencies match with the recent 
experiments. 

PACS numbers: 03.75.Hh, 03.75.Ss, 74.20.Fg 
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Much attention has been focused on Fermion super- 
fluidity realized experimentally using neutral Fermion 
atomic species, ^Li or ^^K 1]. It is achieved by tuning 
the interaction strength via Feshbach resonance. Upon 
sweeping a magnetic field H through the resonance sit- 
uated dX H = 832G in ^Li case, the system exhibits a 
smooth crossover behavior from Bose-Einstein condensa- 
tion at lower side to BCS at high field side. 

Recently two groups [2, |^ have succeeded in produc- 
ing Fermionic superfluid in ^Li with the population im- 
balance where two interacting up and down spin species 
have different particle numbers. In attractive BCS side 
{H > 832G) which we focus on in this paper, the sys- 
tem shows a quantum phase transition as a function of 
the population difference, that is, the relative polariza- 
tion P = \N^ - Ni\/{N^ ^ Ni) Zwierlein et al H 
assigned a critical imbalances P^ = 0.71, beyond which 
normal state becomes stable, in the resonance experi- 
ment dX H c^ 832G, and demonstrated: (i) Pc decreases 
as the temperature T rises, {ii) P^ is proportional to the 
gap value A, that is, Pc ex e~^/^^^l^l ex A (kp Fermi 
wave number), (iii) The spatial profile of the minority 
component exhibits a "bimodal" distribution, which dis- 
appears when the system becomes normal either above 
T > Tc (Tc is the transition temperature) or P > Pc. 
On the other hand. Partridge et al. 3] have found an- 
other transition at P* = 0.09 at the same resonance field 
(832G). 

The problem of the BCS with population imbalance 
posed by these experiments has been addressed in vari- 
ous contexts, ranging from ferromagnetic superconductor 
ErRh4B4 p, heavy Fermions superconductor CeCoIus 
under a field p, to color superconductivity in dense 
quark matter of high energy physics [3 • Among various 
proposals the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
state I3, Q is a prime candidate to describe these exper- 
iments where the sign of the spatially varying gap func- 
tion A(r) changes periodically. This is contrasted with 
the usual BCS state in which A(r) keeps a definite sign 
even though they might spatially vary in some situations. 
Here, we call it the "BCS" state which has a definite sign 



in the gap function. The "FFLO" state changes its sign 
somewhere in the system, even though both states have 
spatially varying gaps due to trapping. 

Prior to the present work, there have been papers de- 
voted to this problem. However, most works are either an 
infinite uniform system without trappi ng 1 1 ( J or does not 
consider possibility of the FFLO state [ullll], except for 
a few J^] . It is crucial to take into account two effects si- 
multaneously because we are considering an intrinsically 
non-uniform finite system. The energy difference of the 
FFLO and the BCS is so subtle because in the FFLO 
solution the sign change occurs only near the surface and 
both spatial profiles are similar at the trap center. 

The purpose of this paper are two- fold: (1) To con- 
struct a generic phase diagram in the plane of T versus 
P and (2) to characterize each phase by the microscopic 
calculation considering a trap. In particular, the origin of 
the observed bimodality in minor component, item (in), 
can be attributed to a characteristic of the FFLO state. 
The derived phase diagram with the Lifshitz point [l^ 
turns out to be quite universal where three second or- 
der lines meet at the tricritical Lifshitz point [l^ llq . 
A new aspect here is that the population in two species 
is the control parameter while in usual condensed mat- 
ter systems the relative shift of the chemical potential 
is controllable. The resulting phase diagram shows that 
at T = the FFLO state is ubiquitous and always the 
ground state for any imbalanced cases up to a critical 
value. It allows us to understand two critical values Pc 
2] and P* 3] at T > in addition to items (i) and (ii). 

It is convenient to use the Bogoliubov-de Gennes 
(BdG) formalism to describe A(r) with population im- 
balance under a trap V{y). Throughout this paper, we 
set fi = /cb = 1. We consider a cylindrical system with 
V(r) = ^Mur'^^ imposing a periodic boundary condi- 
tion with the periodicity bd {d~^ = VMcj) along the 
z-direction. Thus we are treating a three dimensional 
(3D) system, depending on the radius r and z. In the 
current work, the longitudinal trap along the z-axis is 
absent. The BdG equation for the quasi-particle wave 
functions 'Uq(r) and 'i^q(r) labeled by the quantum num- 



ber q is read as follows with the local density of each spin 
state Pa{r) (a =t, |) and attractive interaction g [l7J : 



/Ct(r) A(r) 
A*(r) -/C|(r) 









(1) 



where /Ct,|(r) = -3]^ + V{r) ^ gpi^^{r) - /i|,| with the 
third term being the Hartree term. The self-consistent 
gap equation is given by 

A(r) =gY,u^{r)v*^{r)f{e^) - f A(r)G|;(r). (2) 



Since the chemical potential shift Sja = ja^ — /i| causes 
the braking of the time-reversal symmetry, the sum in 
Eq. (|2|) is done for all the eigenstates with both positive 
and negative eigenenergies [l^- To regularize A(r), we 
have subtracted G^^^ which is the irregular part of the 
single-particle Green's function [lfj |. 

In Fig. 1 we show the phase diagram at T = in the 
plane of P and the energy gap Ao(0) at the trap center 
for equal population case normalized by the Fermi en- 
ergy Ef{0) at the center. It is seen that the critical im- 
balance Pc beyond which the normal state becomes sta- 
ble is linearly proportional to Ao{0)/Ef{0). This shows 
the critical chemical potential difference Spc is given by 
Sfic/Ep = aAo/Ep where the proportional constant a 
delicately depends on the dimensionality of the system 
and/or the Fermi surface shape. For example, one can 
find aiD ^ oo p in one dimension and asr> — 1-41 [2i| 
in 3D Fermi sphere. Since P = |-^ for the normal state 
with 3D Fermi sphere under the assumption Sp <^ Ep^ 
Pc = f<^3D^ — 1-14^, which is also shown in Fig. 1. 
This is changed in the presence of the harmonic trap to 



)<^ho 



Ao(0) 
Ef{0) 



because P 



6fi 
2 EpiO) 



in the normal state. 



Our numerical calculation in Fig. 1 shows aho — 1-26. 

We note that the linear relation Pc ex -^^j^ is observed 
experimentally since the observed phase boundary be- 
tween the normal and superfluid states is approximately 
exponential behavior, namely, Pc oc g-^^/^feFlal ^g^^ p^g ^ 
in Ref. ^). This agreement must be checked further ex- 
perimentally for wider 1//cf|«| region. 

This enhanced critical value Pc given by 6fic/^o = Ci^ho 
far exceeds other known values, such as the so-called 
Pauli-Clogston limiting value S^c/Aq = 1/V^ signaling 
the first order transition from BCS to the normal state, 
or BCS-FFLO unstable point Spc/Ao = 2/7r which cor- 
responds to the one-soliton creation energy 5j. This indi- 
rectly proves that the present solution of the BdG equa- 
tion is stable energetically. 

In order to obtain the observed Pc = 0.7 at the reso- 



Ao(0) 
EpiO) 



0.38, 



nance |2], we can read off from Fig. 1 that 
by performing the naive extrapolation from the weak cou- 
pling limit. This value is compared with the theoretical 
estimate f^ - 0.49 by Bulgac et al HJ. The direct es- 
timation of Pc at unitarity limit is still open to question. 




FIG. 1: (color online) Phase diagram in P vs Ao(0)/^f(0). 
Circles, squares, and triangles show the FFLO, BCS, and nor- 
mal states obtained by our calculation, respectively. Dotted- 
dashed line denotes the result of 3D Fermi sphere case [20|. 
Star is the experimental date on the resonance {H — 832G) 



It is also shown in Fig. 1 that at T = in the stable su- 
perfluid state the order parameter A(r) always exhibits 
the sign change except for equal population (P = 0), 
that is, the FFLO state is stable. The sign of A(r) must 
change to accommodate excess majority species at T = 0, 
which is nothing but FFLO, while at T 7^ the "magne- 
tization" can be accompanied with the non-oscillating 
pairing via thermally excited quasi-particles. The tt- 
shift, i.e., sign change, of the gap function is essence of 
the "topological doping" important for stripes in high 
Tc superconductors ,241. This is also well-known in the 
other physics field [2^; The commensurate (C) charge or 
spin density waves (CDW, SDW) give way to the incom- 
mensurate (IC) ones when adding excess carriers. The 
present problem is precisely analogous to this C-IC prob- 
lem, where IC (C) corresponds to FFLO (BCS). 

We are now in position to characterize the FFLO and 
BCS states. We show two typical examples of these states 
in Fig. 2. In FFLO shown in Fig. 2(a) we can divide 
the density distribution into three distinctive regions (I) , 
(II), and (III) from the center. In "BCS core" region (I) 
the up and down-spin atoms are nearly equal and bal- 
anced. The population of down-spin atoms is enhanced 
and pulled up by mutual attraction. The gap develops 
fully there and the magnetization m(r) = P|(r) — P|(r) 
almost vanishes. In "FFLO" region (II) the gap function 
changes its sign, allowing to accommodate the excess ma- 
jority species. These features give rise to (A) the bimodal 
distribution in the minor component (see a shoulder of 
Pl{r) in lower panel of Fig. 2(a)) and (B) a sharp peak 
structure in m{r) at r/d ^2.7. These features which 
are observed experimentally |3, 1^ come from the sign 
changes of A(r). We note that the ordinary BCS region 
(I) is describable with the local density approximation 
[ill , while the oscillating pairing in (II) results from only 
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FIG. 2: (color online) Spatial distributions of A(r), Pair) 
and m(r) for (a) FFLO state (P = 0.44 and T=0) and (b) 
BCS state (P = 0.22 and T/u;=4.0). Ao(0)/£;f(0)=0.28. cj 
is the trap frequency. 



the full numerical calculation of the BdG equation (1). 
In "complete polarization" region (III) the gap is almost 
vanishing and there is no minority species. Thus the 
complete polarization is attained there. 

These characteristics (A) and (B) are indeed observed 
experimentally |^]. Note that experimental data 2] are 
obtained by the columnar integrated density distribu- 
tions, yet they show prominent bimodality and sharp 
peak structures, implying that the actual three dimen- 
sional features are sharper. 

These characteristics in the FFLO are contrasted with 
those in the BCS shown in Fig. 2(b); The gap has a def- 
inite sign and the two density profiles for two species are 
smooth and scaled to each other. At T > 0, the differ- 
ence /)^(r) — p|(r) = 771 (r) appears by thermal excitations. 
The outer edge contains only the up-spin atoms where the 
gap vanishes. The resulting m(r) has no sharp feature. 
The minority distribution ceases to exhibit the bimodal- 
ity. These features are almost the same as those in the 
normal state given by the Thomas- Fermi profile. These 
features are observed either in the region P > Pc = 0.7 
in Ref. or P < P* = 0.09 in Ref. %,. 

In order to better characterize the FFLO, we display 
series of changes both |A(r)| and vn{r) at T = O in Fig. 3. 
It is seen that as P increases, the periodicity of the os- 
cillations in A(r) relatively stays constant and the num- 
ber of the sign change increases. As for m{r)^ with in- 
creasing P the double peak structure changes into a sin- 
gle peak above Pc, signalling phase transition from the 
FFLO state to the normal state. We can see small and 
faint features in m{r) corresponding to the sign change 
of A(r). 

In Fig. 4 we show the phase diagram in the plane of 
T/Tco versus P/Pco where Tco is the transition tem- 
perature at P = O and Pco is the critical imbalance 
at T = O. This is determined by solving Eqs.(l) and 




FIG. 3: (color online) Spatial distributions of |A(r)| (a) and 
m{r) (b) for various P's. Ao(0)/^f(0)=0.28 and T = 0. The 
normal state appears for P > Pc = 0.57. The dotted-dashed 
lines show the nodes of A(r). 



(2) for various T and P where circle (square) shows 
the FFLO (BCS) state and inverted triangle is the nor- 
mal state. All lines indicate the second phase transi- 
tions which meet at a tricritical point L known as the 
Leung point, or Lifshitz point in more general context 
Q. The BCS-FFLO line starts right from P = T = 0, 
implying that the ground state is always FFLO when 
P 7^ as mentioned above. The BCS appears only 
at a finite T. From Eqs.(l) and (2), it is easy to de- 
rive the equation for Tc, or the boundary between the 
superfiuid and normal state as function of P, namely, 
1 = ^/dr/dr'^*(r)[i^(r,rO - S{t - T')G%{r)]^{r'), 
where 



Kir,r') = J2 



/M-/(£q') 



wq'(r)<,(r')<(r)«q(r').(3) 



The gap function ^(r) ex A(r) is normalized as 
/(ir|^(r)p = 1. Using ^(r) for the BCS state, we es- 
timate Tc. The results are also plotted in Fig. 4 as the 
dashed-dotted line, showing a good agreement with the 
full numerical computation. We have confirmed that Tc 
in FFLO becomes higher than that in BCS beyond the 
Lifshitz point, proving the stability of FFLO over BCS. 
According to the experiment [3 Pc{T) decreases 
as T increases for three magnetic fields {H = 832, 
883, and 924G). More quantitatively, at the resonance, 
Pc{T)/Pco = 0.86 at T/Tf = 0.12 fl, which is indi- 
cated in Fig. 4, estimated within the experimental value 
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FIG. 4: (color online) Phase diagram in T/Tco vs P/Pco 
for FFLO (circle), BCS (square), normal state (triangle) at 
Ao(0)/^f(0) = 0.28. The dashed-dotted line is Tc for BCS 
obtained from Eq. @. Empty circle is the Lifshitz point. Star 
is the experimental data ,4J- 



on the resonance [23: Tco/Tf = 0.27. These agree 
with our results. It is interesting to notice that under 
a fixed temperature, say, T/Tco = 0.3 as P increases, 
BCS changes into FFLO at P/Pco = 0.1 and upon fur- 
ther increasing P, FFLO finally becomes unstable at 
P/Pco — 0.9. It is reasonable that there are two tran- 
sitions observed by Partridge et al. (P* = 0.09) 3] and 
Zwierlein et al. {Pc = 0.71) 2, 4] at the same resonance 
field {H = 832G). The former (latter) is BCS-FFLO 
(FFLO-normal) transition. 

This phase diagram shown in Fig. 4 is quite generic, 
which describes various physical systems, such as CDW, 
SDW, or stripe phase in high Tc superconductors [2J, 
but usually expressed in terms of T vs S/j.. Here the 
population imbalance is a control parameter of the sys- 
tem. This phase diagram in T vs P is qualitatively same 
for other coupling strengths or Ao{0)/Ef{0) or l//ci?|tt|. 
Figure 5 displays a schematic phase diagram in T, P and 
l/kplal of the attractive side only. In T vs I/ZcfI^I plane 
the phase boundary shows Tco ex e~^/^^^l^l. The phase 
boundary in P vs I/ZcfI^I plane is also an exponential 
behavior because Pc ex Aq. 

Since it is rather difficult to distinguish between the 
FFLO and BCS only from the density profiles, we defi- 
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nitely need further experimental probes to identify each 
phase. According to Ref. [3] vortices in the outer region 
(II) and (III) become invisible while only in (I) visible 
vortices are sustained. The quantum depletion of parti- 
cle density at a vortex core only occurs when A/Ep be- 
comes large enough |23,|23- Excess majority species in 
region (II) fill out the vortex core preferentially, making a 
vortex invisible through the density profile measurement 
even though the phase of A(r) winds around. 

In conclusion, we have constructed a generic phase di- 
agram in T vs P and explained various experimental 
aspects, including items (i)-(iii). We have approached 
the strong coupling limit problem on resonance from the 
weak coupling using the BdG formalism by assuming that 
there is no phase transition between them. An under- 
standing of the critical imbalance 0.7 from first principles 
is an outstanding problem, which belongs to future work. 

The authors thank M.W. Zwierlein for valuable discus- 
sions and R. Hulet for useful information. 
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